Supramolecular organization of heptapyrenotide oligomers--an in depth investigation by molecular dynamics simulations.

We present a molecular modeling study based on ab initio and classical molecular dynamics calculations, for the investigation of the tridimensional structure and supramolecular assembly formation of heptapyrenotide oligomers in water solution. Our calculations show that free oligomers self-assemble in helical structures characterized by an inner core formed by π-stacked pyrene units, and external grooves formed by the linker moieties. The coiling of the linkers has high ordering, dominated by hydrogen-bond interactions among the phosphate and amide groups. Our models support a mechanism of longitudinal supramolecular oligomerization based on interstrand pyrene intercalation. Only a minimal number of pyrene units intercalate at one end, favoring formation of very extended longitudinal chains, as also detected by AFM experiment. Our results provide a structural explanation of the mechanism of chirality amplification in 1:1 mixtures of standard heptapyrenotides and modified oligomers with covalently linked deoxycytidine, based on selective molecular recognition and binding of the nucleotide to the groove of the left-wound helix.

Engineering tocopherol selectivity in α-TTP: a combined in vitro/in silico study.

We present a combined in vitro/in silico study to determine the molecular origin of the selectivity of [Formula: see text]-tocopherol transfer protein ([Formula: see text]-TTP) towards [Formula: see text]-tocopherol. Molecular dynamics simulations combined to free energy perturbation calculations predict a binding free energy for [Formula: see text]-tocopherol to [Formula: see text]-TTP 8.26[Formula: see text]2.13 kcal mol[Formula: see text] lower than that of [Formula: see text]-tocopherol. Our calculations show that [Formula: see text]-tocopherol binds to [Formula: see text]-TTP in a significantly distorted geometry as compared to that of the natural ligand. Variations in the hydration of the binding pocket and in the protein structure are found as well. We propose a mutation, A156L, which significantly modifies the selectivity properties of [Formula: see text]-TTP towards the two tocopherols. In particular, our simulations predict that A156L binds preferentially to [Formula: see text]-tocopherol, with striking structural similarities to the wild-type-[Formula: see text]-tocopherol complex. The affinity properties are confirmed by differential scanning fluorimetry as well as in vitro competitive binding assays. Our data indicate that residue A156 is at a critical position for determination of the selectivity of [Formula: see text]-TTP. The engineering of TTP mutants with modulating binding properties can have potential impact at industrial level for easier purification of single tocopherols from vitamin E mixtures coming from natural oils or synthetic processes. Moreover, the identification of a [Formula: see text]-tocopherol selective TTP offers the possibility to challenge the hypotheses for the evolutionary development of a mechanism for [Formula: see text]-tocopherol selection in omnivorous animals.

Molecular origin of piezo- and pyroelectric properties in collagen investigated by molecular dynamics simulations.

Molecular dynamics simulations were used to study the effect of mechanical and thermal stimuli on the electrostatic properties of collagen model helices. Our model sequences were based on glycine proline and hydroxyproline amino acids. We find that longitudinal mechanical strain induces significant variation of the polarization of the collagen fibril. Such a phenomenon is determined by reorientation of the backbone polar groups, which are free to respond to the mechanical solicitation. This non-negligible effect is facilitated by the peculiar folding structure of the collagen helix, which is characterized by the absence of an extended hydrogen-bond network. The stretching/compression of the helix requires a concomitant winding/unwinding motion of the global structure; therefore, the shear components of the stress tensor are the components that most effectively induce structural modification associated to the piezoelectric response. The present calculations also report a pyroelectric response to thermal activation. Model calculations indicate that the pyroelectric effect is dominated by secondary components associated with the piezoelectric tensor.

Common mechanistic features among metallo-beta-lactamases: a computational study of Aeromonas hydrophila CphA enzyme.

Metallo-beta-lactamases (MbetaLs) constitute an increasingly serious clinical threat by giving rise to beta-lactam antibiotic resistance. They accommodate in their catalytic pocket one or two zinc ions, which are responsible for the hydrolysis of beta-lactams. Recent x-ray studies on a member of the mono-zinc B2 MbetaLs, CphA from Aeromonas hydrophila, have paved the way to mechanistic studies of this important subclass, which is selective for carbapenems. Here we have used hybrid quantum mechanical/molecular mechanical methods to investigate the enzymatic hydrolysis by CphA of the antibiotic biapenem. Our calculations describe the entire reaction and point to a new mechanistic description, which is in agreement with the available experimental evidence. Within our proposal, the zinc ion properly orients the antibiotic while directly activating a second catalytic water molecule for the completion of the hydrolytic cycle. This mechanism provides an explanation for a variety of mutagenesis experiments and points to common functional facets across B2 and B1 MbetaLs.

Adaptive protein evolution grants organismal fitness by improving catalysis and flexibility.

Protein evolution is crucial for organismal adaptation and fitness. This process takes place by shaping a given 3-dimensional fold for its particular biochemical function within the metabolic requirements and constraints of the environment. The complex interplay between sequence, structure, functionality, and stability that gives rise to a particular phenotype has limited the identification of traits acquired through evolution. This is further complicated by the fact that mutations are pleiotropic, and interactions between mutations are not always understood. Antibiotic resistance mediated by beta-lactamases represents an evolutionary paradigm in which organismal fitness depends on the catalytic efficiency of a single enzyme. Based on this, we have dissected the structural and mechanistic features acquired by an optimized metallo-beta-lactamase (MbetaL) obtained by directed evolution. We show that antibiotic resistance mediated by this enzyme is driven by 2 mutations with sign epistasis. One mutation stabilizes a catalytically relevant intermediate by fine tuning the position of 1 metal ion; whereas the other acts by augmenting the protein flexibility. We found that enzyme evolution (and the associated antibiotic resistance) occurred at the expense of the protein stability, revealing that MbetaLs have not exhausted their stability threshold. Our results demonstrate that flexibility is an essential trait that can be acquired during evolution on stable protein scaffolds. Directed evolution aided by a thorough characterization of the selected proteins can be successfully used to predict future evolutionary events and design inhibitors with an evolutionary perspective.

Design of HIV-1-PR inhibitors that do not create resistance: blocking the folding of single monomers.

The main problems found in designing drugs are those of optimizing the drug-target interaction and of avoiding the insurgence of resistance. We suggest a scheme for the design of inhibitors that can be used as leads for the development of a drug and that do not face either of these problems, and then apply it to the case of HIV-1-PR. It is based on the knowledge that the folding of single-domain proteins, such as each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES), stabilized by local contacts among hydrophobic, strongly interacting, and highly conserved amino acids that play a central role in the folding process. Because LES have evolved over many generations to recognize and strongly interact with each other so as to make the protein fold fast and avoid aggregation with other proteins, highly specific (and thus little toxic) as well as effective folding-inhibitor molecules suggest themselves: short peptides (or eventually their mimetic molecules) displaying the same amino acid sequence of that of LES (p-LES). Aside from being specific and efficient, these inhibitors are expected not to induce resistance; in fact, mutations in HIV-1-PR that successfully avoid the action of p-LES imply the destabilization of one or more LES and thus should lead to protein denaturation. Making use of Monte Carlo simulations, we first identify the LES of the HIV-1-PR and then show that the corresponding p-LES peptides act as effective inhibitors of the folding of the protease.

Modeling the alpha-helix to beta-hairpin transition mechanism and the formation of oligomeric aggregates of the fibrillogenic peptide Abeta(12-28): insights from all-atom molecular dynamics simulations.

In this paper, all-atom molecular dynamics simulations in explicit solvent are used to investigate the structural and dynamical determinants of the alpha-helical to beta-hairpin conformational transition of the 12-28 fragment from the full length Abeta Alzheimer's peptide. The transition from alpha-helical to beta-structure requires the peptide to populate intermediate beta-bend geometries in which several mainly hydrophobic interactions are partially formed. This is followed by the sudden collapse to ordered beta-hairpin structures and the simultaneous disruption of the hydrophobic side-chain interactions with a consequent increase in solvent exposure. The solvent exposure of hydrophobic side-chains belonging to a sequence of five consecutive residues in the beta-hairpin defines a possible starting point for the onset of the aggregation mechanisms. Several different conformations of model oligomeric (dimeric and tetrameric) aggregates are then investigated. These simulations show that while hydrophobic contacts are important to bring together different monomers with a beta-hairpin like conformation, more specific interactions such as hydrogen-bonding and coulombic interactions, should be considered necessary to provide further stabilization and ordering to the nascent fibrillar aggregates.

Beta-hairpin conformation of fibrillogenic peptides: structure and alpha-beta transition mechanism revealed by molecular dynamics simulations.

Understanding the conformational transitions that trigger the aggregation and amyloidogenesis of otherwise soluble peptides at atomic resolution is of fundamental relevance for the design of effective therapeutic agents against amyloid-related disorders. In the present study the transition from ideal alpha-helical to beta-hairpin conformations is revealed by long timescale molecular dynamics simulations in explicit water solvent, for two well-known amyloidogenic peptides: the H1 peptide from prion protein and the Abeta(12-28) fragment from the Abeta(1-42) peptide responsible for Alzheimer's disease. The simulations highlight the unfolding of alpha-helices, followed by the formation of bent conformations and a final convergence to ordered in register beta-hairpin conformations. The beta-hairpins observed, despite different sequences, exhibit a common dynamic behavior and the presence of a peculiar pattern of the hydrophobic side-chains, in particular in the region of the turns. These observations hint at a possible common aggregation mechanism for the onset of different amyloid diseases and a common mechanism in the transition to the beta-hairpin structures. Furthermore the simulations presented herein evidence the stabilization of the alpha-helical conformations induced by the presence of an organic fluorinated cosolvent. The results of MD simulation in 2,2,2-trifluoroethanol (TFE)/water mixture provide further evidence that the peptide coating effect of TFE molecules is responsible for the stabilization of the soluble helical conformation.

Understanding the determinants of stability and folding of small globular proteins from their energetics.

The results of minimal model calculations indicate that the stability and the kinetic accessibility of the native state of small globular proteins are controlled by few "hot" sites. By means of molecular dynamics simulations around the native conformation, which describe the protein and the surrounding solvent at the all-atom level, an accurate and compact energetic map of the native state of the protein is generated. This map is further simplified by means of an eigenvalue decomposition. The components of the eigenvector associated with the lowest eigenvalue indicate which hot sites are likely to be responsible for the stability and for the rapid folding of the protein. The comparison of the results of the model with the findings of mutagenesis experiments performed for four small proteins show that the eigenvalue decomposition method is able to identify between 60% and 80% of these (hot) sites.